Bogomol'nyi, Prasad and Sommerfield configurations in smectics 
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It is typical in smectic liquid crystals to describe elastic deformations with a linear theory when 
the elastic strain is small. In smectics, certain essential nonlinearities arise from the requirement of 
rotational invariance. By employing the Bogomol'nyi, Prasad and Sommerfield (BPS) decomposi- 
tion and relying on boundary conditions and geometric invariants, we have found a large class of 
exact solutions. We introduce an approximation for the deformation profile far from a spherical in- 
clusion and find an enhanced attractive interaction at long distances due to the nonlinear elasticity, 
confirmed by numerical minimization. 
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Exact solutions to nonlinear model problems offer deep 
insight into systems not at first apparent from simple 
linear approximations. Liquid crystals with smectic-A 
order, consisting of a stack of regularly spaced fluid lay- 
ers, have particularly interesting nonlinear elastic prop- 
erties that arise from considerations of rotational invari- 
ance. The combination of these nonlinearities and ther- 
mal fluctuations is the source of anomalous elasticity in 
smectics Q . More recently, there has been someprogress 
in the nonlinear theory of screw dislocations |2| and in 
the structure of isolated edge dislocations Q. Nonethe- 
less, the consequence of these essential nonlinearities has 
been mostly ignored in the study of deformations induced 
by defects such as edge dislocations and inclusions em- 
bedded between the layers 0, 0, 0| • In this letter we 
show that a large class of defect configurations can be 
probed by appealing to their implied boundary condi- 
tions and geometrical invariants. We find, in particular, 
that though the layer displacements (shown in Fig. 1) 
can greatly differ between the linear and nonlinear the- 
ories, that the interactions between edge dislocations is 
unchanged. 

Recently, Brener and Marchenko have developed a so- 
lution for the deformation around a single edge disloca- 
tion that takes into account some of these nonlineari- 
ties 0, and that differs significantly from the deforma- 
tion profile derived from linear elasticity, even very far 
from the defect where the elastic strain and layer cur- 
vature is small. This deformation profile has been con- 
firmed in an experiment by Lavrentovich and Ishikawa 
in a cholesteric finger texture. Despite this success, it is 
difficult to understand how to use the solution to gain in- 
sight on the effects of nonlinearities on more complicated 
deformations (for example, curved or multiple edge dis- 
locations). 

We have developed a large class of one-dimensional ex- 
act solutions which include the solution of Q as a limit- 
ing case. This larger class includes, among others, certain 
multiple edge dislocation configurations. The method of 
solution takes advantage of the fact that the free energy 



FIG. 1: The layer displacement of four equally spaced, 
parallel edge dislocations in the nonlinear (solid) and lin- 
ear (dashed) elasticity theories. The Burgers vector for all 
dislocations is 20A and the layers shown are at z = X, 
z — 4A, 8A, • • • . Note that the linear and nonlinear solutions 
differ by as much as 4A and that the nonlinear solution relaxes 
more rapidly. 



density for any one-dimensional deformation can be writ- 
ten as the sum of a perfect square and a total deriva- 
tive. This form also arises in the study of vortices in 
superconductors and kinks in nonlinear field theories . 
Expressing the free energy in this form reveals that the 
total deformation energy is always bounded from below 
(the Bogomol'nyi bound) and that so-called BPS config- 
urations saturate this bound. These BPS configurations 
satisfy a nonlinear differential equation of reduced order, 
which can be transformed into a linear equation through 
the Hopf-Cole transformation. Thus, we are able to cal- 
culate the layer deformation and energy of these edge 
dislocation configurations exactly and ponder the effects 
of superposition. 

This same technique produces approximate solutions 
to deformations that are not one-dimensional and that, 
in many cases, have lower energy than the deformation 
profile derived from linear elasticity. In particular, we 
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study the deformation profile far from an embedded in- 
clusion in a lamellar phase, and the asymptotic inter- 
action energy of two such inclusions. These results are 
relevant, for instance, to the stud y of focal conic defects 
with positive Gaussian curvature [HJ, UJ$ ■ 

The free energy for a smectic-A liquid crystal is 



F = — 



d 3 3 
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+ A 2 (V-n) 
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where B is the bulk modulus, A is a length defined by 
the fact that X 2 B is the bending modulus of the smec- 
tic layers, n = (z — Vw) /|z — Vu| is the unit normal to 
the smectic surface, and w(x) is the Eulerian displace- 
ment field. The first term accounts for the compression 
strain between layers while the second term arises from 
the bending stiffness of the layers. The latter is man- 
ifestly rotationally invariant; the former can be shown 
to be rotationally invariant as well. This nonlinearity is 
required and symmetry dictates the relative numerical 
factor of i between the first- and second-degree terms. 

In the limit of small strains, it is typical to retain 
only terms quadratic in derivatives of u. This leads 
to a linear elasticity which misses much of the essen- 
tial physics. Instead, we make the approximation that 



d 



(V ± u) <C 1, where the notation Vj^ has been 



used to represent the two dimensional gradient in the 
plane of the layers. In the following, we will write diU to 
represent the in-plane derivatives and continue to write 
d z u for the derivative along the layers at infinity. In this 
limit, we rewrite the free energy density as: 
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We shall see that the energy density is the sum of a per- 
fect square and a total derivative under the appropriate 
conditions. This sort of form for the free energy has 
been used successfully to study the fluctuations of two- 
dimensional smectics 0] and to study shape changes in 
vesicles 0]. We use it here to study defects. 

We can simplify this expression further by considering 
another total derivative, the Gaussian curvature K = 
\V-[n (V • n) — (n • V) n], the divergence of a vector 
To lowest order in the displacement field u, K K — 
\di (diud 2 u — djudjdiu). We decompose the vector into 
longitudinal and transverse components: 



diudj u — djudjdiu 



-eijdjcj) - diip 



(3) 



for scalar fields <j) and ip. Using this expression, we are 
able to rewrite the last terms in J5J as: 

d z uV\u = di (d z udiu) - i<9 2 (diuf (4) 

3(V ± u) 2 V 2 ± u = -AKu + diidiUdjudju) 

—2di(ueijdjcf) + udiip) (5) 



We see immediately that if K — (as would be the case 
for a one-dimensional deformation such as an edge dislo- 
cation), the free energy is a perfect square plus a series 
of surface terms, i.e. the BPS form. Note that since 
J d 2 xK is a constant depending on the boundary con- 
ditions (similar to the topological invariant J d 2 xK), a 
constant shift of u merely shifts the zero of energy so 
translation invariance, u — > u + e, is preserved. 

A result that follows immediately from this observa- 
tion is that the energy of any deformation with K = 
is bounded from below by the boundary contributions 
to the free energy. This bound is saturated by certain 
minima of F, namely those for which the perfect square 
term vanishes and which, consequently, satisfy a first- 
order differential equation in contrast with the second- 
order equation which would follow from extremizing the 
free energy. If we consider a deformation that depends 
only on z and x then K — and so if 



x - (d x u) 2 T Xd 2 u = 0, 



(6) 



then the free energy will depend only on the boundary 
terms. Through the Hopf-Cole transformation S± — 
e ±u/2\ we gee fasti S± satisfies the diffusion equation 
d z S± = Xd 2 S±. Solving equation 10 we find that the 
positive branch of the theory describes deformations for 
z > and the negative branch describes deformations 
for z < 0. By "gluing" the two solutions together, we 
can determine the free energy merely by considering the 
boundary at z = 0. 

For a single edge dislocation, the boundary conditions 
are S± — > 1 as x — > — oo and S± — > e ±h / 4A asu +oo, 
where b is the Burgers vector of the dislocation. We find 



S+ = A+F, 



2VXz 



(7) 



where E(a;) = (vr)" 1 / 2 f*^ dtcxp(-t 2 ), C± = 1 and 

A± = (e b / 4A — l) • It is straightforward to invert the 
Hopf-Cole transformation to find agreement with the 
edge dislocation deformation calculated in Q. Since 
boundary terms are essential to our analysis, we note that 
this solution can also be obtained through the boundary 
condition, u±(x,z = 0) = ±bO(x)/2 where Q(x) is the 
step function. 

The contribution to the line tension from the defor- 
mation of the layers results from a single boundary term 
from each region, z > and z < 0. This can be cal- 
culated by approximating the step function at z = 
with the integral of a Gaussian of width £, resulting in 
the expression t = F/L y = BXb 2 /3y/irt;. This has the 
same form as the line tension of an edge dislocation in 
the quadratic approximation for the free energy up to a 
rescaling of the microscopic length £. The reason for this 
is simple: were we to rewrite the energy of the linearized 
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theory in the same form as J3J), we would find only one 
surface term. However, in either the linear or nonlinear 
theories the boundary is the xy-plane and the free energy 
comes entirely from total derivatives with respect to z. 
Thus only the (common) second term in 0} contributes 
to the energy and this BPS configuration of the nonlinear 
theory has precisely the same energy as the same solution 
of the linear theory, though, as seen in experiment , the 
actual layer deformations are quite different. Though our 
solutions are absolute minima for our imposed topology, 
we have not explored the possibility of modifying the 
defect geometry while maintaining the same boundary 
conditions at infinity. For instance, it is possible that fo- 
cal conic defects may introduce additional boundaries in 
the bulk which could lower the energy further. A more 
complete analysis will be necessary to address this issue 

Because S satisfies a linear equation, we can superpose 
solutions to calculate the interaction between two, par- 
allel edge dislocations. For two parallel edge dislocations 
at z = and positions X\ < xi and with Burgers vectors 
b\ and hi respectively, we have 



sgn(z)u(x,z)/2A 
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The layer deformation for a set of four, equally spaced 
parallel edge dislocations located at z — is shown in 
Figure^ 11 comparison to the deformations derived using 
a linear elastic theory. Since only one boundary condition 
can be specified (at z = 0), one must explicitely check 
that the BPS solution satisfies the remaining boundary 
condition, d z u — > as z — > oo in this case. Due to 
this, multiple edge dislocations solutions in different lay- 
ers are not BPS configurations unless they are infinitely 
far apart. 

The energy of a multiple edge dislocations can be de- 
termined in an analogous manner to that of a single edge 
dislocation. Since the free energy is entirely determined 
by the same boundary term in both cases, the free en- 
ergy of many edge dislocations must be equal to the free 
energy of those edge dislocations in the harmonic approx- 
imation. As a result, the interaction energy between edge 
dislocations in the same layer must also agree. This result 
is in stark contrast to the interactions of screw disloca- 
tions. There, in the linear theory, screw defects interact 
exponentially, while in the nonlinear theory they interact 
via a power-law repulsion [l^ . 

Up to this point, solving Eq. (0 led to exact min- 
ima of the free energy. We could interpret that equation 
as equating the strain, u zz — d z u — h (V ± u) 2 , to the 
bending AH = AV^m which is physically compelling: it 
equates the two terms in the energy in an attempt to 
minimize the frustration between the two ;18|. For more 




FIG. 2: The layer displacement (points) at r = for the 
boundary condition of u(r,z = 0) = 9A exp(— r 2 /25A 2 ), com- 
paring the numerical minimization with the nonlinear (solid 
line) and linear (dashed line) elasticity. In the inset is a com- 
parison of the peak shape at z = 50A for all three cases. Ex- 
cept for a region about r = 0, the BPS approximation tends 
to be closer to the numerical minimization. 



general boundary conditions (K ^ 0), the free energy 
is not in the BPS form. Nonetheless, if the additional 
term, —ABXuK, is small compared to the energy scale 
set by the curvature, B(W\u) 2 (which is the only other 
energy scale), we can presumably treat that term as a 
perturbation. In order to test this hypothesis, we have 
compared the layer deformation due to a Gaussian in- 
clusion with the results from a numerical minimization 
using a conjugate-gradient algorithm (results shown in 
Figure EJ). The layer displacement at r = falls between 
the BPS solution and the harmonic solution. Though the 
harmonic solution is more accurate for small r, overall 
the BPS approximant is closer to the numerical solution. 
This is not surprising: because Ku 3> 1 near r = we 
should not expect the BPS solution to be valid there. 
As we choose boundary conditions with smaller values of 
K, the BPS solutions become more and more accurate. 
Further study is needed to understand the precise role of 
the uK in the failure of the BPS configurations to mini- 
mize the energy |15|. It is often the case that "near-BPS" 
solutions are remarkably good approximants fl6| and it 
appears to be true here as well. 

We consider a spherical inclusion of radius R located 
at z = which implies the boundary conditions u± (r, z = 
0) = ±V-R 2 — r 2 Q(R — r) for z > and z < respec- 
tively, where r — \J 'x 2 + y 2 , and u±(r, z) vanishes for 
r —>■ oo and z — * ±oo. Solving u zz — AH, we have 



u(x, z) — 2Asgn(z)< 1 + 
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where x are coordinates in the plane of the layers. 

Far from the inclusion the strain in small and we per- 
form a multipole expansion for r, z 3> R and find that 
the monopole term dominates. We have 

u(r,z) = 2Aln(\ + ^e- r2/4A ^ , (10) 

the deformation of a monopole in the region z > 0, where 
/3 = R (e R / 2X - H) + 2A (1 - e fl / 2A ). Due to its expo- 
nential dependence on the radius of the inclusion, we 
see that for R > 2A, (3 ~ Rexp(R/2X) becomes very 
large and so (|10(1 differs greatly from the linear solution 
U{ = f3e~ r / 4Xz /z. Higher order moments may also be 
calculated explicitly for a spherical particle. 

With the success of the single monopole solution, 
we turn to the interaction of two inclusions that are 
well separated. We expect the deformation to be well- 
approximated by a superposition of BPS monopole con- 
figurations. For two monopoles of size fti and /?2 sepa- 
rated by a distance 2xo and both located at z = 0, the 
deformation profile can be found by applying superpo- 
sition. By expanding in powers of Pi, we can estimate 
the interaction energy. The lowest order term is, in prin- 
ciple, quadratic in (3i and identical to the interaction in 
the harmonic approximation and higher order terms arise 
from the nonlinearity of the elasticity: 

t I g V n J 

Xq 

20TTl3 1 f3 2 (512/3? + 245/3i/? 2 + 512/%) 
+ Slag + " ' ' 

Notice that the term quadratic in the monopole moments 
vanishes and thus, as in the case of screw dislocations 
[T?! the interaction energy is essentially nonlinear. This 
expression for the energy compares reasonably with our 
numerical estimates. 

We have introduced a method to calculate deforma- 
tions of a smectic liquid crystal in a way that takes into 
account nonlinearities in the elasticity, as well as their en- 
ergies. While the method is asymptotically exact for the 
one-dimensional deformations, we find that it produces 
deformations which have lower free energies than those 
of the usual quadratic approximation. This allows us to 
calculate the asymptotic interaction energy between two 
large smectic inclusions. In future work we will study the 



implications of nonlinearities of dislocations with both 
screw and edge components, and the nonlinear elasticity 
of the full smectic free energy. 
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